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We report an experimental study of the primary bifurcation in electrically-driven convection in 
a freely suspended film. A weakly conducting, submicron thick smectic liquid crystal film was 
supported by concentric circular electrodes. It electroconvected when a sufficiently large voltage V 
was applied between its inner and outer edges. The film could sustain rapid flows and yet remain 
strictly two-dimensional. By rotation of the inner electrode, a circular Couette shear could be 
independently imposed. The control parameters were a dimensionless number TZ, analogous to the 
Rayleigh number, which is oc and the Reynolds number TZe of the azimuthal shear flow. The 
geometrical and material properties of the fllm were characterized by the radius ratio a, and a 
dimensionless number V, analogous to the Prandtl number. Using measurements of current-voltage 
characteristics of a large number of films, we examined the onset of electroconvection over a broad 
range of a, V and TZe. We compared this data quantitatively to the results of linear stability theory. 
This could be done with essentially no adjustable parameters. The current-voltage data above 
onset were then used to infer the amplitude of electroconvection in the weakly nonlinear regime by 
fitting them to a steady-state amplitude equation of the Landau form. We show how the primary 
bifurcation can be tuned between supercritical and subcritical by changing a and TZe. 



I. INTRODUCTION 

Phenomological amplitude equation models are impor- 
tant tools in the study of symmetry-breaking, pattern- 
forming bifurcations. They have proven particularly 
useful at the weakly nonlinear level. While their broad 
application stems from the common symmetries that are 
shared by many different systems, it is the coefficients 
or specific parameters in amplitude equations that dis- 
tinguish each system. It is an interesting question to ex- 
amine how these coefficients, especially those of the non- 
linear terms, change as the properties of the system are 
varied. There are only a few experimental systems, such 
as Rayleigh-Benard convection, Taylor vortex flow, and 
electrohydrodynamic convection in nematic liquid crys- 
tals, for which these coefficients have been quantitatively 
determined. 

Three-dimensional, nonlinear systems are prone to 
develop complicated spatial and temporal patterns 
even when only weakly nonequilibrium. j^] The spatio- 
temporal patterns are often the result of successive sym- 
metry breaking bifurcations. It is useful to study pattern 
formation in low-dimensional systems which are close to 
equilibrium but have little symmetry so that there are 
only a very limited set of symmetry-breaking bifurca- 
tions available. In general, one seeks the most complex 
dynamics that can be realized in as simple and restricted 
a system as possible. In our system, annular electro- 
convection, we exploit the strict two-dimensionality of a 
submicron smectic A liquid crystal film. The lower di- 
mensionality greatly reduces the variety of possible pat- 
tern states and so makes it easier to experimentally study 
the rich nonlinear properties of the basic pattern. 



In this paper, we report an experimental study of the 
properties of the bifurcations in annular electroconvec- 
tion with a variable circular Couette shear. Our primary 
focus is to investigate the variation of the cubic nonlinear 
coefficient of an amplitude equation as the experimental 
parameters are systematically changed. This coefficient 
determines the saturation behaviour of the convective 
flow velocity. We begin by establishing that the exper- 
imental system is well-described by a theoretical model 
that reduces to an amplitude equation near the onset of 
the pattern-forming instability. We test the adequacy of 
the theoretical model by direct comparisons between the 
experimental results and theoretical predictions at the 
linear level. |^ We then proceed to describe how the co- 
efficient of the lowest order nonlinear term depends on 
the experimental parameters. Interestingly, we find that 
the coefficient can be made to change sign and thus that 
the primary bifurcation can be tuned between subcritical 
and supercritical. In this paper, we show from symme- 
try considerations that the amplitude equation that de- 
scribes the sheared case is the complex Ginzberg-Landau 
equation. Elsewhere, we have derived this result 

directly from the full electrohydrodynamic equations us- 
ing a multiple-scales analysis. [|| Calculation of the pa- 
rameter dependence of the coefficients of this equation is 
beyond the scope of this paper, but in future, we expect 
to be able to compare the experimental results presented 
here with theory at the weakly nonlinear level. 

Our system consists of a freely suspended annular 
smectic A film which can be driven out of equilibrium 
by electrical forces and can be independently subjected 
to a shear flow as shown in Fig. It was previously 
described in Rcfs. S and 0. The film was suspended 
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between concentric circular electrodes and was driven by 
a voltage difference between the electrodes. The inner 
electrode could be rotated about its axis, thereby im- 
posing a circular Couette shear. The flow pattern that 
emerges when the film is driven sufficiently hard is re- 
ferred to as electroconvection and consists of an array of 
counter-rotating vortices. Our primary probe of the elec- 
troconvection amplitude was the excess electrical current 
carried convectively by the flow. The annular geometry 
and electrical nature of the experiment make it possible 
to study radially forced convection under steady shear, a 
combination more often associated with geophysical flows 
than with small laboratory systems. On the other hand, 
the accurately two-dimensional nature of the flow makes 
it far simpler than most planetary analogs, and thus more 
amenable to theoretical investigation. 

The driving force originates from the interaction of the 
radial electrical field and the surface charge density that 
develops on the film's free surfaces. Below the onset of 
convection, the film's electrical state is independent of the 
imposed shear in the film. It is thus straightforward to 
superpose shear flows with the radial driving, as these are 
characterized by independent control parameters. The 
simplest such shear flow is steady azimuthal or circular 
Couette flow. The annular geometry is naturally periodic 
so that the shear flow is closed and leads to a net mean 
flow around the annulus. 

The addition of a shear flow to the annular electrocon- 
vection alters the symmetries of the base state of elec- 
troconvection. When shear is absent, the base state is 
invariant under azimuthal rotation and reflection in any 
vertical plane containing the rotation axis through the 
center of the annulus. The electroconvective state is then 
stationary and appears with the spontaneous breaking of 
the azimuthal invariance. Q When sheared, the reflection 
symmetry of the base state is not present. When driven, 
the pattern once again breaks the azimuthal symmetry, 
but since the reflection symmetry is absent due to shear, 
the pattern is free to travel azimuthally in the direction of 
the mean flow. |^] In addition, as we show in detail below, 
the shear flow alters the primary bifurcation, making it 
hysteretic. 

The equilibrium properties of freely suspended liquid 
crystal films have been extensively studied; see Ref. Q for 
a recent review. Smectic liquid crystals consist of layers 
of orientationally ordered long molecules which readily 
form suspended films. In smectic A, the average orien- 
tation of the long axis of the molecules, and hence the 
optic axis, is normal to the layer plane. The layers are 
of uniform thickness and within each layer the distribu- 
tion of molecules is isotropic. Smectic A exhibits two- 
dimensional isotropic fiuid properties in the layer plane 
while fiows perpendicular to the layers are strongly inhib- 
ited. Other material properties such as the electrical con- 
ductivity and dielectric permitivity are also isotropic in 
the plane of the layers. Uniform suspended smectic films 
are always an integer number of smectic layers thick and 
while they flow readily, they seldom change thickness; 
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FIG. 1. The geometry of the annular film. Circular Cou- 
ette shear is produced by rotating the electrode holding the 
inner edge of the film. 

they are robustly two-dimensional. Unlike soap solutions 
P,plo| , smectic liquid crystals have very small vapor pres- 
sures, hence the film can be enclosed in an evacuated 
environment. The reduction in ambient pressure leads 
to a proportional reduction in the air drag to which the 
film is subjected. This system thus has several at- 
tractive features for the study of bifurcations in simple 
nonlinear systems that may be described by amplitude 
equation models. Our study greatly extends and com- 
plements previous work on electroconvection in isotropic 

smectic A 10. 



liquids |11[-|17[, nematics 
and higher smectic phases 




The paper is organized as follows. In Section II A 



we 

describe the results of a linear stability analysis of our 
system. In the process we define the various dimension- 
less para meters that characterize the experiment. In Sec- 
tion II B, we use symmetry arguments to justify the form 
of the relevant amplitude equation. Our description of 
the experimental apparatus and protocol follows in See- 
In Section IV A we outline our data analysis 
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procedure. We present experimental results that per- 
tain to linear analysis in Section IV B. An experimental 



determination of the cubic coefficient of the amplitude 
equation is found in Sections [V C and IV D . Section ^ 
discusses the relationships between our results and those 
for other similar systems and section VI presents a brief 
summary and conclusion. 



II. THEORETICAL BACKGROUND 

In this section, we briefiy review the linear and weakly 
nonlinear theory that is relevant to our data analysis and 
results. 

A linear stability analysis of annular electroconvection 
with circular Couette fiow was reported in Ref. The 
theory is constructed with a number of simplifying as- 
sumptions. 1^,0 The fluid fllm is treated as an annular 
sheet of inner radius r^, outer radius Tq and thickness s. 
A schematic is shown in Fig. |l|. An important geometric 
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parameter is the radius ratio a = ri/vo- The film width 
d = To — Ti is assumed to be much greater than s. The 
viscous fluid is assumed to flow only in two-dimensions, 
and be incompressible, Ohmic and of uniform electrical 
conductivity. We denote the fluid density by p, its molec- 
ular viscosity by 77 and its electrical conductivity by a. 
Only the charges due to free surfaces are included and 
bulk dielectric effects inside the film are neglected. 
The electrodes are assumed to be of negligible thickness, 
and fill the rest of the plane not occupied by the annular 
film. Outside the plane of the film, we assume an empty 
space of dielectric permittivity eq- Couette shear in the 
theoretical model is imposed by specifying the appropri- 
ate velocity boundary conditions on the edges of the film. 
Linear theory predicts the position of the onset of con- 
vection and the degree to which onset is suppressed by 
the shear. These are discussed in detail in Ref. Pl and 



are further analyzed in Section IV B below 



From very general symmetry considerations, we can 
deduce the form of an amplitude equation which is valid 
in the weakly nonlinear regime just above onset. ||l|,||j4| 
Here, a complex-valued amplitude describes the slowly- 
varying magnitude and phase of the spatially periodic 
physical fields, for example the fiuid velocity. We find 
an expression of the well-known Landau form. The 
coefficients of the amplitude equation can in principle 
be calculated from the underlying elcctrohydrodynamic 
equations; a calculation of this kind is in progress and 
will be reported separately The magnitude of the 
complex amplitude is directly related to the quantity we 
measured experimentally, the total current carried by the 
film. We can thus fit current-voltage data to the real part 
of the amplitude equation to experimen tally deter mine 
its coefficients. This is done in Sections IV C and [VP 
below. 



A. Outline of linear stability theory 

The theoretical treatment of Ref. |^ examined the lin- 
ear stability of the two-dimensional annular fluid fllni 
when it is subjected to a voltage V at the inner elec- 
trode of the annulus while the outer electrode is held 
at ground potential. In addition, the theory allowed for 
the rotation of the inner edge of the annulus at angular 
frequency w, while the outer edge is held fixed. The ro- 
tation drives a Couette shear flow in the film below the 
onset of convection. It can be shown that any arbitrary 
rotation of the inner and outer edges can be transformed 
to a rotation of only the inner edge by moving to the 
appropriate rotating coordinates. [p|j30|] Such a coordi- 
nate transformation has no effect on the dynamics as 
long as the flow is strictly two-dimensional. The theoret- 
ical model is completely specifled by the radius ratio a 
and three additional dimensionless parameters: 



7^ = -2-^ 
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and -Re = i—— . (2.1) 



The control parameter 7?. is a measure of the external 
driving force which is proportional to the square of the 
applied voltage V while P is a ratio of the time scales 
of electrical and viscous dissipation processes in the fllm. 
The Reynolds number TZe is a measure of the strength of 
the applied shear, which is regarded as a second control 
parameter. It has been established that the instabil- 
ity leads to a one-dimensional pattern of m vortex pairs 
where m is the azimuthal mode number. Linear stabil- 
ity predicts the value of the critical control parameter 
TZc and the critical mode number rric at which the fllm 
becomes marginally unstable. TZc and rric are in general 
functions of a, V, and TZe. A detailed discussion of the 
calculated values of TZc and rric under various conditions 
is given in Ref. When TZe = {i.e. with no applied 
shear), it was found that TZc and ttlc are independent of 
T'. In the following, these critical values for zero shear 
will be denoted 7^° and m°. 

When TZe > 0, it was found that TZc > TZc', for any a 
and "P, i.e. that the shear always suppresses the onset of 
convection. It is convenient to measure the relative 
degree of suppression for a given a and 7-" in terms of the 
reduced quantity 
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which is a monotonically increasing function of TZe for 
all "P. In Section IV B below, we extend the discussion 
in Ref |^ by comparing these results to data for several 
different a, and for a wide range of T'. 



B. Nonlinear theory: The amplitude equation 

The amplitude equation appropriate to our system fol- 
lows from symmetry considerations. We defer a 
detailed calculation of its numerical coefficients to future 
work. Q 

In the absence of shear, the base state solution of 
the elcctrohydrodynamic equations is invariant under az- 
imuthal rotations and refiections through any plane per- 
pendicular to the annular plane and containing its cen- 
ter. Just above onset, we assume that all the physical 
fields are nonaxisymmetric and proportional to a rapidly- 
varying spatial oscillation e""^ , where m is an azimuthal 
mode number. This requires that the differential equa- 
tion for the slowly-varying complex amplitude Am of 
mode m, be unchanged under the transformations 
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(2.3) 



(2.4) 



In the above, 9 is an azimuthal angle, and the overbar 
denotes complex conjugation. The most general ampli- 
tude 
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FIG. 2. Schematic of the electrical and mechanical parts 
of the apparatus, showing the annular electrodes and film 
drawing assembly in side view. 



equation that is invariant under these symmetry opera- 
tions has the Landau form, 

TdtA^ = eA,n - g\A,n\^A.m + /i|A„/A„ - ... , (2.5) 

where r, g, and h are real-valued coefficients. The small 
parameter e is a reduced control parameter given by 



Tic 



-1, 



(2.6) 



and is a measure of the distance from threshold. The 
amplitude equation is accurate for small e and describes 
a bifurcation from the Am = state (e < 0) to the Am ^ 
state with 2m vortices (e > 0). A sharp bifurcation 
occurs at e = 0. 

If the fluid is subjected to a circular Couette shear, the 
base state is only invariant under azimuthal rotations, 
Eqn. 2.3. In this case, the general amplitude equation 



takes the complex Landau form 

T{dt - iaim)Am = e(l + ica)Am 
+ h{l + ic:i)\A„ 



- g{l + iC2)\Am\'^ Am 

\^Am ~ ... , (2.7) 



where ajm is the imaginary part of the eigenvalue of the 
unstable mode m at onset and the coefficients cq, C2, C3 
are real-valued. Eqn. 2.7, which here is simply deduced 
from symmetry considerations, has been rigorously de- 
rived for annular electroconvection with shear from the 
basic electrohydrodynamic equations. j|] In general, it 
describes a Hopf bifurcation to a pattern with a complex 
amplitude Am which travels in one azimuthal direction. 
In both Eqns. 



2.5 and 2.7 



we have omitted terms 
involving the azimuthal gradients of Am, which if in- 
cluded would give the so-called Ginzburg-Landau equa- 
tion. ^ While such terms are potentially quite interest- 
ing, they do not appear to be directly relevant to inter- 
preting our current-voltage data. Our annular system 



is azimuthally periodic and thus lacks lateral boundaries 
at which Am — > where gradients would be important. 

Thus, any gradients in Am would have to arise from 
purely dynamical effects, such as Eckhaus or other phase 
instabilities 

From current-voltage data, we can extract the reduced 
Nusselt number n, which is a global dimensionless quan- 
tity defined as the ratio of current carried by convection 
to that by conduction. The reduced Nusselt number n is 
related to the electric Nusselt number M hy n = M — 1. 
One can show (2^J^ that the amplitude can be scaled so 
that n = l^mp. 

To solve for the real and imaginary parts of the am- 
plitude equation, let Am{t) = A(t)e'*^^ where A{t) and 



<^{t) are real. Substitution into Eqn. p. 4 gives 



TdtA ^(A- gJ^ - hA^ 
T(dt^ - aim) = (Co - gc2A^ + ... 



(2.8) 
(2.9) 



We will show in Section [V A how the raw current-voltage 
data can be transformed into measurements of e and n, 
and hence be used to measure the real amplitude A via 
a nonlinear fit procedure. We focus on extracting the 
coefficient g of the cubic nonlinearity. The sign of g de- 
termines whether the pitchfork bifurcation to electrocon- 
vection is supercritical ( "forward" , g > ) or subcritical 
( "backward" , g < ). A tricritical bifurcation occurs 
when g — 0. The results for g are discussed in Sec- 
tions |ivq and |IV4 



III. EXPERIMENTAL APPARATUS 

In this Section we describe the liquid crystal sample 
and the experimental apparatus. We performed a large 
number of precise current-voltage measurements on elec- 
troconvecting annular films with a wide range of radius 
ratios and applied shears. In overall structure, the ap- 
paratus consisted of the annular electrodes which were 
housed in a vacuum chamber, an electrometer circuit 
and a stepper motor, both under computer control. The 
chamber also allowed films to be drawn under vacuum 
and inspected under reflected white light to check thick- 
ness uniformity. A schematic is shown in Fig. 0. 



A. The liquid crystal 



As in previous experiments [pyqj20|-|24| , our study used 
smectic A octylcyanobiphenyl (8CB). 8CB has a smec- 
tic A phase between 21°C and 33.5°C. The electrical 
conductivity of "pure" 8CB is due to several ionic im- 
purities of varied and unknown concentrations. In order 
to control the nature of the ionic species contributing to 
the electrical conductivity, we doped the 8CB with tetra- 
cyanoquinodimethane (TCNQ), a material believed to 
form charge transfer complexes with the host, so that the 
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dominant impurity species was the dopant. To prepare 
the doped material, TCNQ was dissolved in acetonitrile 
and added to the 8CB sample. The acetonitrile was then 
evaporated in a vacuum oven while warming the mixture 
so that the 8CB was in its isotropic phase. The sam- 
ples used had concentrations of 2.96 x 10^^, 1.11 x 10^"' 
and 7.62 x 10"^, of TCNQ by weight. Experiments with 
significantly higher or lower dopant concentrations were 
found to be less useful due to irreproducible non-ohmic 
behaviour below the onset of convection. [B2| 



B. the annulus 

The annular electrodes were constructed out of stain- 
less steel. The inner electrode was a circular disk of di- 
ameter 2ri. The outer electrode was a circular plate of 
diameter 9.00 cm with a central hole of diameter 2ro. 
The outer electrode was 0.73 ± 0.01 mm thick. By us- 
ing several pairs of inner and outer electrodes, we con- 
ducted experiments at six different radius ratios between 
a = 0.33 and a — 0.80. We used inner electrodes with 
radii ranging between — 3.60 ± 0.01 mm and — 
5.26 ± 0.01 mm. The radii of the outer electrodes were 
between Tq = 5.57 ± 0.01 mm and Tq = 11.25 ± 0.01 mm. 

The annulus was housed in a vacuum chamber which 
was evacuated by a rotary pump. The air was evacuated 
slowly so as to prevent vigorous air flows that may cause 
the film to rupture. The surrounding air was pumped 
down to an ambient pressure between 0.1 — 5.0 torr. At 
these pressures, the mean free paths of iV2 and O2 are 
in the range 0.5 — 0.05 mm. This was comparable to the 
film width d, so that the air drag on the film was negli- 
gible. 1^ All experiments were performed at the ambient 
room temperature of 23 ± 1°C, well below the smectic 
A-nematic transition at 33.5°C for undoped 8CB. 

The inner electrode was made to rotate about its axis 
by means of a high precision stepper motor operating at 
25600 steps per revolution. The motor was located out- 
side the vacuum chamber and was connected through a 
rotating seal. The inner electrode was adjusted to rotate 
true to its axis to within 50^m at angular frequencies up 
to a; = 67r rad/s. 

C. Drawing the film and determining its thickness 

A motorized film-drawing assembly was housed in the 
vacuum chamber. The film was spread across the annular 
gap between the electrodes using a stainless steel razor 
blade inclined at 25°. The blade was moved away from 
the annulus while voltage was applied, to avoid perturb- 
ing the electric fields near the film. 

The thickness of the film was determined from the in- 
terference color of the film under reflected white light 
using a low power videomicroscope. Using standard col- 
orimetric functions |33-p5t| a color-thickness chart was 



calculated for 8CB. Since smectic films are constrained 
to be an integer number of smectic layers thick, the film 
color can be used to identify the film thickness measured 
in smectic layers. Each layer of smectic A 8CB is 3.16 nm 
thick. | |3^ Most of the experiments were performed with 
films between 25 and 85 layers thick. Over most of the 
middle of this range, the film thickness can be deter- 
mined to within ±2 layers, while close to the ends of the 
range a more conservative estimate of ±5 layers was used. 
The colorimetric method worked well for film thicknesses 
where the strong interference color can be unambiguously 
matched to a color chart. Very thin films that appeared 
black and thicker, nearly white films were not used. 

During the course of an experiment, the film color was 
monitored to check that it remained uniform in thick- 
ness to within ±1 layer. Films drawn with non- uniform 
thicknesses, left to themselves, tend to anneal to become 
uniform in thickness and hence uniform in color. The an- 
nealing process can be accelerated by electroconvecting 
and shearing the films. 

Current-voltage runs during which the film thickness 
spontaneously became non-uniform were abandoned. A 
small number of visualization experiments were per- 
formed on films which had a nonuniformity in thickness of 
about ±2 smectic layers, or about 5% of s. Usually these 
nonuniform films had two thicknesses and in refiection 
displayed two nearby interference colors. The advection 
of the thickness nonuniformities was used to visualize the 
flow pattern and thus to determine the azimuthal mode 
number m by counting vortices. For simplicity, some of 
these observations were done at atmospheric pressure. 
This visualization method, while somewhat crude, was 
preferrable to using suspended smoke particles, which 
tend to perturb the conductivity of the film. [p3| 

D. Current- voltage measurements 

Since the current transported through the film was pi- 
coamperes in magnitude, particular care had to be ex- 
ercised to avoid stray currents. The inner rotating elec- 
trode was electrically isolated from all other components 
by a highly insulating sleeve. The current was carried 
onto the rotating electrode by a set of silver-graphite slip 
rings inside the vacuum chamber. A Keithley model 6517 
electrometer was used both as a voltage source and a pi- 
coammeter. The 'high' of the programmable dc voltage 
source of the electrometer was connected to the rotating 
inner electrode, which is the only part of the entire appa- 
ratus that is not at ground potential. The outer annulus 
was connected to the input of the electrometer which 
is held at very nearly ground potential. The outer an- 
nulus was electrically isolated from true ground by tefion 
washers which served to eliminate leakage currents which 
would otherwise be added to the signal. The rest of the 
apparatus was grounded. Electrical noise was reduced by 
shielding the electrodes and most of the 
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FIG. 3. Representative current-voltage characteristics for 
radius ratio a = 0.467 in the absence of shear (a) and when 
strongly sheared (b). Note the different scales. 

experimental components in a large Faraday cage which 
doubled as the vacuum chamber. Low noise triaxial ca- 
bles and feedthroughs were used to connect the outer 
electrode to the electrometer. 

An example of a current- voltage characteristic in a film 
without shear is shown in Fig. It consists of data 
obtained for incremental and decremental voltages. The 
current-voltage characteristic clearly shows two regions: 
one for voltages smaller than a critical voltage Vc and 
one for voltages greater than Vc- The critical voltage that 
separates these two regions is in the vicinity of the kink in 
the current-voltage characteristic. In the regime V < Vc, 
the current is linearly dependent on the voltage and the 
film is ohmic. Experiments in 8CB with significantly 
higher and lower concentrations of TCNQ show, at least 



initially, non-ohmic current-voltage characteristics even 
for V <Vc, due to electrochemical effects. 

We estimated Vc before each run then defined an ap- 
proximate reduced control parameter e — (T^/Vc*^**)^ — 1. 
Current-voltage data was then obtained by making vari- 
able steps in voltage in such a way as to make equal in- 
crements in e. We used several e step sizes, which became 
finer near threshold. Similarly, we used a variable waiting 
time after each step, which allowed for longer relaxation 
times closer to Vc- Very long relaxation times were not 
feasible due to the drift of the electrical conductivity. |3^] 
After the wait time, between 100 — 200 measurements of 
the current, each separated by 25 milliseconds, were av- 
eraged. The error in the average was taken to be the 
standard deviation of the mean or 1%, the reading error 
of the picoammeter, whichever was greater. Each run 
consisited of incrementing e up to a predetermined max- 
imum and then decrementing it to zero voltage again. 
Once the data had been acquired, more precise e values 
at each voltage were found using drift-corrected values of 
Vc from a fit procedure described in the next section. 

When the film was sheared by rotation of the inner 
electrode, it was allowed at least 30 seconds after a 
change of shear rate to attain a steady state before the 
current-voltage characteristics were obtained. In all runs, 
the shear rate was established and then held fixed during 
a current-voltage sweep. Fig. ^ displays a representa- 
tive current-voltage characteristic in a film under shear. 
The principle effects of the shear are to suppress the on- 
set of convection and, eventually, to make the primary 
bifurcation hysteretic. 



IV. EXPERIMENTAL RESULTS 

With the exception of some simple flow visualization 
which used slightly nonuniform films, all our results were 
obtained by fitting current-voltage data taken with uni- 
form films. It is possible to fit the data in such way as to 
extract all the unknown material parameters in a nearly 
model-independent way, as well as to correct for small 
drifts in conductivity. 

The first task was to determine the dimensionless pa- 
rameters relevant to each current-voltage sweep: the ra- 
dius ratio a, the dimensionless ratio P and the Reynolds 
number TZe. We can then properly non-dimensionalize 
the current-voltage data and fit it to an amplitude equa- 
tion model. The details of this fit are described in the 
section 



IV A below. 



In Section [V B , we consider the features of the data 



that are predicted by the linear stability theory devel- 
oped in Ref. For zero shear, these are the critical 
mode number m[?, and the critical voltage V^. We use 
the latter to fix the last unknown material parameter, 
the viscosity rj. We also compare our data to the linear 
stability prediction for the shear-suppression of the onset 
of convection. This was done for six values of a, and for 
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wide ranges of V and TZe. 

In the weakly nonlinear regime, the value of the coeffi- 
cient of the cubic nonlinearity, g in the a mplit ude equa- 
tion, was of primary interest. In Section IV C we report 
results for g at various a and V for TZe — 0, that is, in 
the absence of shear. We find that g can become slightly 
negative for small a, so that the bifurcation is weakly 
backward. The result for g in the limit TZe = 0, a — » 1 is 
compared to the theoretical result for rectangular films 
p5| . Finally, in Section IV D| we describe the results for g 
in sheared films, where TZe > 0. We find that increasing 
TZe has a strong effect on the nature of the bifurcation, 
driving g negative and making the bifurcation backward. 



A. Fitting tlie current-voltage data 

Except for a small drift, the film is ohmic below the 
onset of convection. We used fits in this range to de- 
termine some important material parameters in order to 
scale the data. 

Each voltage-current measurement in the ohmic 
regime, (V, /), constitutes an experimental determination 
of the film's conductance, c = I/V. For a film of radius 
ratio a, thickness s and conductivity a, the conductance 
is given by 



2TTas 
ln(l/a) 



(4.1) 



Interestingly, the conductance is independent of the size 
of the film, i.e. independent of or Tq. Since a is 
merely a geometrical parameter, measurements of c are 
effectively measurements of as, eliminating the need to 

determine thes e se parately. 

Using Eqn. 4.1, T', defined in Eqn. 2.1, can be ex- 
pressed in terms of c as 



V = 



eo?7 



27reo?7 



1 



pasd p(ro - Ti) In (l/a) 



(4.2) 



where d = Tq — r^. Similarly, the Reynolds number TZe — 
pujrid/r] can be calculated from the measured angular 
frequency lu of the inner electrode in rad/s, given p and 
T]. The density p of 8CB at room temperature |Q is 
1.0 X lO'^ kg/m'^. In both T' and TZe , the only remaining 
undetermined material parameter is the viscosity r). We 
fixed this parameter using a combination of ex perim ental 
and theoretical results, as described in Section IV B . The 
result was 77 = 0.18 ± 0.03 kg/ms. 

Thus, for each current-voltage sweep, we were able 
to determine 7^ and TZe. The drift in conductivity | |3^ ] 
caused a drift in c and hence in P over the course of a 
sweep. Furthermore, c is only known directly during the 
ohmic parts at the beginning and end of each sweep. We 
associated a mean T' with each sweep by averaging over 
c data before and after. This results in ~ 10% errors in 
■p for any one sweep, and a slow, uncontrolled evolution 
of T' from sweep to sweep. 



The whole range of current-v oltag e data is fit to the 
amplitude equation model. Eqn. 2/7 gives the full, time- 
dependent complex coefficient amplitude equation re- 
quired by sy mmetry. Only the real and steady state part 
of Eqn. 2.7, given by 



eA - gA^ - hA' + f ^ , 



(4.3) 



is required to model the current-voltage data. Here 
we h ave truncated at the quintic order and augmented 
Eqn. 2.7 with a phenomenological field term /. This term 
models the rounding of the bifurcation due to nonideal 
systematic effects which are slightly symmetry-breaking, 
resulting in an "imperfect" bifurcation. ||l[ For all of the 
fits, we found / ^ 1. 
As discussed in section 



II B 



the amplitude A which ap- 
pears in Eqn. 4.3 is related to the reduced Nusselt number 
n by 



Icond 



1 



1 = 



(4.4) 



while the reduced control parameter e is given by 

.„^-l.(XV-l. (4.5) 



TZc 



Vc 



Eqns 4.3-4.5 contain five parameters, c, Vc, 5, /i, and 
/. While c has a nearly constant initial slope and 
is marked by a relatively obvious kink in the raw (/, V) 
data, the amplitude A is i ndi rectl y de duced via the pair 
of transformations Eqns. 4.5 and 4.4 which are nonlin- 
ear. This amplitude is in turn fit using Eqn. [4.3| , which 
is again nonlinear. Thus, the parameters (g, /i, /) are 
rather distantly related to the raw (/, V) data. If g < 0, 
the bifurcation is hysteretic, and A is multivalued over 
some ranges of e. Also, the nature of the model neces- 
sarily involves several fit parameters which are not inde- 
pendent. Consequently, the determination of these pa- 
rameters is muc h m ore di fficult than V^. As discussed 
in Sections IV C| and IVD , {g,h, f) can be influenced by 
small systematic effects in the data leading to scatter 
which is larger than the statistical uncertainties in the 
fits. Nevertheless, the general trends are robust. 

The drift of the film's conductivity, and hence its con- 
ductance c, introduces some additional complications. 
Recall that the onset of convection occurs when the con- 
trol parameter TZ equals a critical value TZc given by 



TZc 



Vc - -^/^z, 



(4.6) 



Since the electrical conductivity is not constant, the value 
of Vc required to calculate e slowly changes during the 
course of an experiment. Combining Eqns. 4.1 and 4.6 
yields 



Vc(c) 



' s?]TZc ln(l/a) c 



27reo2 



(4.7) 
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FIG. 4. The amplitude A and the reduced Nusselt number 
n = vs. the control parameter e for a film with radius 
ratio a = 0.64. The solid and dashed lines are a fit to the 
Landau amplitude equation. 
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FIG. 5. The amplitude A and the reduced Nusselt number 
n = A^ vs. the control parameter e for radius ratio a = 0.80 
with an applied shear. The solid and dashed lines are a fit to 
the Landau amplitude equation. 



We corrected for the drift in c by tracking its value during 
the ohmic parts of the sweep and using a Hnear interpo- 
lation during the portion of the sweep where the film is 
convecting. In this way, each voltage and current mea- 
surement was non-dimensionalized with its own value of 
c and Vc{c) to produce fittable values of e and A. We 
also transformed the errors in the current measurements 
A/ into amplitude errors AA. 

We approached the fitting problem in the following 
bootstrap fashion. It was easy to ascertain bounds on 
Vc by inspecting the current-voltage characteristics, see 
for example Figs. ||a and b. Each current-voltage char- 
acteristic was scrutinized and two voltage intervals were 
chosen. The first interval contained the critical voltage at 
which the conducting state became unstable to the con- 
vecting state. The second interval contained the voltage 
at which the convecting state became marginally unsta- 
ble to the conduction state. Guesses of Vc in both these 
intervals were chosen at random using a uniform deviate 
random number generator ||39| ] . The conductances c were 
then determined for the ohmic regimes and from a linear 
interpolation for the convecting portion, as described ear- 
lier. The raw current-voltage data w as t hen transformed 
into e and A ± AA and fit to Eqn. 4.3 by varying only 
the three parameters h and / in a weighted Levenberg- 
Marquardt nonlinear fitting procedure which minimized 
X^. We then did a Monte-Carlo optimization of the ran- 
domly chosen Vc- [^,0. The best Vc was taken to be 
the one that minimized over the Monte Carlo sample. 
The uncertainty in Vc was the standard deviation of the 
uniform deviate on the constraint interval. Correspond- 
ing to this best Vc were the three best fit parameters 



g, h and /. The uncertainties in g, h and / were esti- 
mated by a Monte Carlo data decimation step with Vc 
constrained to its best fit value. 

For several reasons, it was necessary to restrict the fit 
to the neighborhood of e ~ 0. In the first place, am- 
plitude equation models are only rigorously valid in the 
limit e ^ 1, although in practice they have been found 
to apply over a more extended range. Secondly, re- 
stricting the range of e reduced the impact of the resid- 
ual, uncorrected component of the conductivity drift. In 
many cases, the drift effects were sufficiently large that 
only the data acquired with increasing e were fit. 

Figure |4| shows a typical result of transforming current- 
voltage data to reduced Nusselt numbers n and ampli- 
tudes A and fitting to Eqn. 4.3. The transition from 



conduction to convection is continuous and the best fit 
parameter g > 0, indicating that the bifurcation is super- 
critical. Figure || shows the analogous result for a film 
under a strong shear. For this case, the transition from 
conduction to convection is discontinuous and g < 0; the 
bifurcation is subcritical. In all fits, we found h > 0, and 
< / ^ 1. We performed a large number of such fits 
and surveyed the dependence of g on a, P, and TZe. In 
addition, the fits provided determinations of the critical 
voltage Vc as a function of a, V, and TZe via the bootstrap 
process descri bed above . These results are discussed in 
Sections llV C| and llVDl below. 
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B. Tests of linear stability theory 

This section compares the experimental measurements 
with the predictions of the hnear stabihty theory given 
in Ref. §]. 

The simplest feature of linear theory that can be com- 
pared with experiment concerns the critical mode num- 
ber for zero shear, mj?. As mentioned in Section [II, some 
experiments were performed in films with slight thickness 
nonuniformity. This permitted the qualitative visualiza- 
tion of the flow field and a quantitative determination of 
the mode number m of the stationary pattern of vortices. 
Unless the bifurcation is strongly subcritical, we expect 
m = mj? close to onset. The observed zero shear mode 
number was in excellent agreement with predictions of 
linear stability analysis. Table | summarizes the results. 
Rapid rotation of the vortex pattern and larger hysteresis 
prevented a systematic study of m under shear. Quali- 
tative observations confirm the general prediction that 
mciTle > 0) < m°, i.e. that shear reduces the number of 
vortices. 

The primary theoretical result of Ref. is the pre- 
diction of the critical voltage Vc required for the onset 
of electroconvection. Vc is given for the general case by 



Eqn. 4.6. Even though the viscosity 77 is not well-known, 



we can test the zero shear theory for various thicknesses 
s and conductivities a using the following scheme. De- 
noting the zero shear value of Vc by VJ?, we write 



arjs 

17" 



(4.8) 



Using Eqn. Eqn. [4 
nicntly as 



can be expressed more conve- 



An^eo'aiVc')^ia) 
= 77c 



(ln(l/a))27^0( 



a 



(4.9) 



Written in this way, Eqn. 4.9 expresses a proportionality 
between and in which the viscosity rj is the only un- 
known parameter. Consistency with this proportionality 
over a wide range of parameters serves as a test of the lin- 
ear theory for Vc as well as a determination of rj. There is 
one caveat: this analysis is not entirely experimental but 
requires the theoretical value of TZc{a). The quantity v, 
which we refer to as the scaled critical voltage, was found 
as follows. For each set of cu rrent voltage data, the fit 
procedure outlined in Section [VA was used to deduce 
a critical voltage Uj* and conductance c at onset. The 
film thickness s was deduced from the color of the film. 
Using the measured radius ratio a and Eqn. [4.1| , a was 
calculated. Finally, the numerical result for 7?.^ (a) was 
used to find which was plotted against c^. Figure ^ 
shows the results obtained from 228 current- voltage runs 
at six different a, numerous different s and conductivi- 
ties in the range 5.9 x 10"^ < cr < 8.4 x 10^'^fi^^m"^ 
Consequently, the range of V is very broad. 
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FIG. 6. The scaled critical voltage vs. the square of the 
conductance for various films at six different a. The solid line 
is a one-parameter fit with a single constant of proportionality 
which determines the viscosity 77 = 0.18 ± 0.03 kg/ms. 



Within some scatter. Fig. ^ exhibits the predicted pro- 
portionality over several decades. Hence, the theory 
properly accounts for the scaling of the critical voltage 
with respect to the film thickness and radius ratio. 

A single-parameter linear fit to v'^ = rjc^ gives ij — 
0.18 ± 0.03 kg/ms. This is a reasonable value for the 
viscosity; while it has not been independently measured, 
it is believed to be of order of 0.1 kg/ms. ^] 

We now turn to the case of nonzero shear. The main 
prediction of the linear stability analysis is that the onset 
of electroconvection is suppressed by Couette shear. This 
suppression is both a and V dependent. Returning to 
Eqn. 2.2, the degree of suppression e is given by 



nc{a,ne,v) 

n{a) 



- 1 = 



Vc{a,ne) 
V^{a) 



- 1 



(4.10) 



We used the two equivalent expressions for e in Eqn. 4.10| 
to calculate the suppression theoretically and experimen- 
tally. Vc{a,TZe) was found from the nonlinear fit, along 
with the conductance c for that particular sweep. It was 
necessary to correct Vc{a) for the drift of c in order to 
calculate i. This was done using values of Vc{a) taken 
from zero shear sweeps performed before and after each 
sheared sweep. The variation of Vj?(a) with c was mod- 
elled with a linear fit which was used to find the drift- 
corrected value of Vc (a) . In this way each experimental 
value of e could be be computed for a constant c. The 
uncertainty in e was dominated by the uncertainty 
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FIG. 7. A comparison between the measured sup- 
pression e vs. TZe and the predictions of local 
and nonlocal theory based on Ref. [3]. In (a) 
a = 0.47. The different symbols denote the P-quartiles: 
13.3 < o < 15.4 < • < 17.5 < □ < 19.6 < A < 21.7. 
The theoretical lines 

are for the mean V = 16.3 of the data. Similarly, in (b) 
a = 0.64, 29.1 < o < 37.1 < • < 45.2 < □ < 53.2 < A < 61.2 
and the theoretical lines are for the mean V = 45.2. Note the 
very different scales. 



cal curves are calculated for the mean value of the range 
of V spanned by the data. For a discussion of the dis- 
tinction between the local and nonlocal approximations, 
see Ref. g. 

Note that the ranges of TZe for these two a are differ- 
ent by a factor of ten and that the suppressions are also 
very different. It is quite interesting that a rather mild 
shear {TZe ^ 3) suppresses the onset of electroconvection 
so strongly that e ^ 14 {i.e. TZc = 15TZ^). We have also 
studied the suppression at a = 0.33,0.56,0.60, and 0.80. 
The results are similar to Fig. |^. The suppression data 
covers the range 10 < T' < 130 and < TZe < 3. In 
each case the agreement between theory and experiment 
is fair, with the theory slightly underestimating the sup- 
pression in most cases. The upper bound on TZe could 
easily be extended using increased rotation rates, while 
smaller T' would require larger, thicker films. 

The degree of agreement shown in Figs. |^ is essentially 
independent of the value of rj. Recall that rj was deter- 
mined by a single parameter fit to Eqn. |4.9| . Since the 
J] dependence in the TZe scaling of both the theory and 
the experiment are proportional to 1 /ri, any change in rj 
multiplies both by the same factor. H This simply re- 
sults in a rescaling of the TZe axis in Figs. |^a and b, with 
no change in the quality of the agreement. 

We close this section with some brief speculation on 
what might account for the scatter in Fig. ^ and for the 
remaining discrepancies between theory and experiment 
in the suppression curves. 

The most likely sources of these discrepancies are the 
imperfect geometry of the film and its finite thickness. 
The 2D theoretical model assumed that the flow velocity 
is independent of position over the thickness of the film. 
This may be inexact for thicker films. Since the electrical 
forcing is localized near the free surfaces, it seems likely 
that the surfaces are preferentially driven in thick films so 
that the motion is not accurately 2D. The film's edges are 
also imperfect because small wetting layers unavoidably 
exist on the circumferences of the inner and outer elec- 
trodes. These may produce electrical or velocity bound- 
ary conditions that are not exactly those assumed by the 
theory. In general, however, the linear stability theory 
works remarkably well considering its rather simple as- 
sumptions. 

C. Coefficients of the Cubic Nonlinearity without 
Shear 



in the drift correction to . The conductance c for each 
e can be used to calculate its associated value of T'. Thus, 
we arrive at sets of experimental values of e{a, TZe, V) for 
fixed a and TZe and ranges of "P. These can be com- 
pared in detail to the theoretical results for suppression 
as a function of these parameters given in Ref. 1^,^. 
Figures |^a and b show this comparison at a — 0.47 
and a = 0.64 respectively. In each case the theoreti- 



In this section, we consider the experimental results in 
the nonlinear regime, beginning with the unsheared TZe = 
case. These are expressed in terms of the cubic Landau 
coefficient g of the amplitude equation, Eqn. 4.3, Even 



with TZe fixed at zero, we will show that this coefficient 
is an interesting and nontrivial function of the remaining 
dimcnsionless quantities, T' and a. 

The radius ratio a is a fixed geometric parameter for 
each run, determined only by the choice 
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FIG. 8. (a) Measurements of the coefficient of the cubic 
nonlinearity g for zero shear vs. V for (a) a = 0.33 and (b) 
a = 0.64. 



of the electrode dimensions. The V appropriate to each 
film was deduced after each run as described in Sec- 
tion [IV A| . Since V is proportional to the conductivity, 
which exhibits a slow drift, a wide range of V could be 
investigated. While they are independent in principle, a 
and V are not easy to separate exp erim entally. This prac- 
tical constraint derives from Eqn. O in which V oc d"^, 
where the width of the film d — Tq — ri. In practice, d is 
largest at small a, so the runs with the smallest V occur 
for small a. 

The nonlinear coefficient g{a^'P) was determined as 
described in Section [V A . Figures ||a and b show g as 

The scatter that is 



a function of V for two different a. 
manifest in these plots exceeds the statistical uncertainty 



of the fit. As discussed in Section IV A, the scatter orig- 
inates from systematic effects due to the non-ideal fea- 
tures of the experiment. Nevertheless, the gross features 
of the V dependence of g can still be extracted. 

At a = 0.33, the measurements explored the range 
2 < V < These were the smallest V reached in the 
experiment. It is clear from Fig. ^a that g increases with 
7^, and passes through zero for small enough V. Thus, 
for a = 0.33, we find that the bifurcation from the con- 
duction state to the m = 4 vortex state is subcritical 
{g < 0) for P ~ 5 and supercritical > 0) for 'P ^ 5. 
Near V = 5 we pass a tricritical point. It is interesting 
to observe that in each case the bifurcation involves the 
same two symmetry states, yet its subcriticality depends 
on V. 

For all the other, larger a investigated, g was found to 
be independent of P. This is illustrated in Fig. ^ for a = 
0.64 and was also true for a = 0.47, 0.56, 0.60, and 0.80. 
However, for all of these cases, P was never less than 10 
and can be as large as several hundred. Thus, it is unclear 
whether the subcriticality at a = 0.33 is due to the small- 
ness of P or the smallness of a, or some combination. It 
would be interesting to determine whether g also becomes 
negative for large a asP decreases. The regime of large a 
and small P, while accessible in principle, would require 
significantly larger electrodes or much thicker films. All 
electroconvection experiments on freely suspended films 
in the rectangul ar g eometry have reported supercritical 
bifurcations [p0|-p3[. However, these experiments were 
also at large P and so the possibility of a subcritical bi- 
furcation in rectangular films cannot be excluded. 

The case of large P is easier to understand. In the 
governing electrohydrodynamic equations (see Ref. [^), 
the parameter analogous to the Prandtl number only ap- 
pears as the inverse, P~^, multiplying certain nonlinear 
terms. Hence, it is perhaps not surprizing that g becomes 
independent of P for P ^ 10, where these terms become 
negligible. 

In order to examine how g depends on a, we removed 
the P dependence by averaging over data at various P. 
For the five largest a, g is independent of P over broad 
ranges oi P. A weighted average of g was obtained for 
these cases. For a = 0.33, where some P dependence was 
evident, only g values in the narrow range 2.1 < P < 4.4 
were averaged. 

In this averaging the systematic scatter in g was 
treated as random. It is thus likely that the true un- 
certainty in the average values of g will be much larger 
than the standard deviation of the mean. The ultimate 
test of this procedure lies in the comparison of the av- 
eraged values of g with theoretical predictions. As we 
describe below, the limited comparison we are able to 
make at present is not unfavourable. The line in Fig. ^3 
shows the average value of g for the plotted data. All of 
the results for g{a) are tabulated in Table || and plotted 
in Fig. |. 
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FIG. 9. Measurements of the coefficient of the cubic non- 
linearity g for zero shear vs. the radius ratio a. The dashed 
line shows the theoretical value for g in the limits V ~* oo 
and a ^ 1 from Ref. [25]. 
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FIG. 10. Measurements of the hysteresis width 5e vs. V 
for radius ratio a — 0.33. 



It is clear from Fig. ^ that, overall, g{a) increases with 
a. As discussed above, the larger a data are also associ- 
ated with larger V . At present, no theoretical predicion is 
available for g at arbitrary a and V . However, we expect 
g to approach a limiting value quite rapidly as a — > 1. 
This limit corresponds to an unbounded lateral geome- 
try in which the film is a strip of fluid suspended at its 
long parallel edges by two semi- infinite plate electrodes. 
In this limit the discrete azimuthal mode numbers m are 
replaced by a continuous wavenumber k = m/f, where 
f = {ri + T^jl is the mean radius. From linear theory, 
we found that the critical parameters 7?.^ and fcj? = m^/f 
at a = 0.80 are already very close to the limiting values 
for a ^ 1. @ 

This limiting behaviour is reasonable given the large 
dimensionless circumference for large a. The relevant as- 
pect ratio A is the circumference at the mean film radius, 
7r(ri + To), divided by the width of the film d = Tq — r;, 
so that A = 7r(l a)/(l - a). At a = 0.80, A ~ 28. 
Almost all the experiments performed in the rectangular 

geometry had A - 10 ||^-||| and were well-modelled by 
the theory for an unbounded strip for which A = oo. It is 
thus resonable to expect 5 at a = 0.80 to be close to its 
limiting value for a = \. Similarly, as the large a data 
also involve very large values of we may also employ 
the theory in the V ^ 00 limit. 

Weakly nonlinear analysis of unsheared electroconvec- 
tion in the 'plate' electrode geometry was presented in 
Ref. for the case V = 00. The result of that analysis, 
g = 2.842 is shown in Fig. || It is in good agreement 
with a reasonable extrapolation of the data to the a — > 1 
limit. 



When g < 0, the onset of convection becomes hys- 
teretic and the quintic term with coefficient h in Eqn. 4.3 
becomes significant. j4|] Since < / ^ 1, the width of 
the hysteresis loop Se in e is given by Se = g^ /4:h. When 
(7 > 0, (5e = 0. The hysteresis width with zero applied 
shear is plotted in Fig. |l^ for a = 0.33 and various V. 
Note that the hysteresis vanishes for V ^ 5 and even 
when non-zero, is always rather small, Se ^ 0.02. This is 
in contrast with the sheared case, discussed in the next 
Section, for which 5e ^ 0.1. 



D. Coefficients of the cubic nonlinearity with shear 

In the presence of shear, the coefficient g of the cubic 
nonlinearity is strongly dependent on the Reynolds num- 
ber TZe of the Couette flow imposed in the base state. 
We find that the shear drives g negative, producing a 
hysteretic Hopf bifurcation to convection in the form of 
travelling vortices. In this Section, we describe the TZe 
dependence of g, which, as in the previous Section, is also 
a function of the radius ratio a and the dimensionless ra- 
tio P. 

Figure ^ shows a representative example of the TZe 
dependence of g, in this case for a = 0.47. At TZe — 0, 



the value of g is found, as described in Section IV C , by 
averaging over a range of V, which was always > 10. In 
the case shown in Fig. 0, g{Tle = 0) = 1.64 ± 0.06 with 
a mean V of 16.3. The plotted g data for TZe > had 
13.3 <V < 21.7. 

Our main result is that TZe plays the part of a second 
control parameter which allows the primary bifurcation 
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FIG. 11. Measurements of the coefficient of the cubic non- 
hnearity g vs. the Reynolds number TZe at a = 0.47. The 
different symbols denote data in various quartiles of V. 



FIG. 12. The hysteresis width Se vs. the Reynolds number 
TZe for radius ratio a = 0.47. The symbols denote the same 
quartiles of V as in Fig. 11. 



to electroconvection to be tuned between super- and sub- 
critical. For a = 0.47, the bifurcation is supercritical at 
TZe = and weakens as the Reynolds number increases. 
The bifurcation becomes tricritical at TZe ^ 0.2 and is 
subcritical thereafter. The subcriticality at first deepens 
with increasing Reynolds number, until at TZe ~ 0.85 a 
minimum value of g — —3.7 is reached. For TZe > 0.85, 
the bifurcation remains subcritical but g becomes an in- 
creasing function of the Reynolds number. For the range 
of Reynolds numbers investigated the bifurcation does 
not become supercritical again. 

There is some scatter in the data, nonetheless, the 
overall trends are clear. The systematic deviations are 
comparable to those in Figs. ||a and b. The results were 
qualitatively similar for a — 0.56, 0.60, 0.64, and 0.80. 
The tricritical Reynolds number at which g — will be 
denoted TZer. The value of TZer varys for different a 
and V, but the general features of the g vs. TZe curve 
are preserved over a wide range of parameters. Table III 
lists the values of TZex under various conditions. 

The variation of TZex is probably due to the combi- 
nation of changing both a and T', rather than to the 
variation either parameter separately. These are difficult 
to experimentally disentangle and the parameter space is 
large. For zero shear, g was found to be independent of T' 

for "P ~ 10, but under shear there is insufficient data to 
draw many conclusions about the "P dependence of g. If 
linear theory is any guide, we expect there will in general 
be a greater T' dependence in the sheared case. It was 
established in Ref. that the linear theory result for TZc 
was independent of "P for TZe = 0, 



but weakly dependent on P for TZe ^ 0. It appears from 
the tabulated results that as P and a decrease, TZer in- 
creases. 

The minimum value of g as a function of TZe also vary 
for different a and P. We will refer to the coordinates 
of this special point as gmin and TZemin- Table IV lists 
the minimum values observed. A look at Fig. |ll| shows 
that the minimum value is only known up to the density 
of data and to within the scatter in the whole plot. The 
minimum is a local one and is observed within ranges of 
TZe which had different maxima for different a. 

Since the parameter space for the data is defined by 
several parameters it becomes difficult to make meaning- 
ful comparisons of the experimental results when more 
than one parameter changes. However, one fortunate 
comparison can be gleaned from Table At a = 0.47, 
the film had P = 15.3 while at a = 0.80, P = 12.0. Since 
the radius ratios are very different and the P are not, 
it is not unreasonable to directly compare the values of 
gmin and TZemin for these cases. It is evident that the 
bifurcation is much more strongly subcritical at a = 0.80 
than at a ~ 0.47, and the values of TZemin are also very 
different. 

The coefficient of the quintic nonlinearity h as function 
of TZe is studied in Ref. |Q. Figure [ij shows the the 
hysteresis width 6e = g^/4:h as a function of TZe for a = 
0.47. The hysteresis when TZe > is much larger than 
that described in the previous Section which was found 
for TZe = at small a and P. 

The Nusselt number of the convection is independent 
of the Hopf frequency {i.e. the travelling rate of the vor- 
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tex pattern). It would be interesting to measure this 
frequency, which is governed by the imaginar y pa rt of 
the complex Ginzberg-Landau equation, Eqn. 2.9, and 
depends on e, a, TZe and V. This would however require 
visualization of the flow pattern or some spatially local- 
ized probe of the amplitude. 



V. DISCUSSION 

Because of its rather simple symmetry and forcing 
scheme, annular electroconvection under applied shear 
has many similarities to other systems and models. In 
this section, we discuss these similar systems with a view 
to putting annular electroconvection into a general con- 
text. 

Electroconvection is obviously analogous to buoyancy- 
driven thermal convection. In fact, it can be shown that 
the linear stability problem for surface charge driven elec- 
troconvection in thin films reduces to the 2D Rayleigh- 
Benard problem if the nonlocal coupling between fields 
and charges is neglected. Detailed calculations show 
that this "local" approximation is an accurate one, even 
when applied in the weakly nonlinear regime. 

There have been some theoretical studies of 3D 
Rayleigh-Benard convection (RBC) in the presence of 
a plane Couette shear flow Q. The theory assumed 
the usual theoretical geometry of a fluid layer confined 
between infinite, perfectly conducting horizontal planes. 
Our annular system approaches the plane one in the 
a — > 1 limit. (The a < 1 annular geometry has finite 
extent and curvature; the implications of these are dis- 
cussed below.) Linear stability analysis of a plane Cou- 
ette base state to RBC reveals stability differences be- 
tween transverse roll disturbances (with axes perpendic- 
ular to the shear flow), and longitudinal roll disturbances 
(with axes parallel to the shear flow). Longitudinal-roll 
disturbances have identical stability properties to un- 
sheared RBC, and are always more unstable than the 
transverse-roll disturbances. In fact, longitudinal-roll 
disturbances have stability properties that are indepen- 
dent of any uni-directional shear flow along the axis of the 
rolls. In our 2D system, these more unstable longitudinal 
roll modes do not exist and the vortices we see correspond 
to truly 2D transverse rolls. Thus, we observe the analog 
of a state that would normally be preempted by longitu- 
dinal rolls if the geometry were not constrained. 

According to linear theory, transverse-roll disturbances 
in unbounded, sheared RBC, like our vortices, exhibit 
suppression, or added stability due to the shear, under 
plane Poiseuille or plane Couette or any mixture of these 
two flows. The onset Rayleigh number for transverse 
rolls is a monotonically increasing function of the shear 
Reynolds number, similar to what we found for electro- 
convection. Furthermore, the critical wavenumber of the 
most unstable transverse disturbance was found to be a 
monotonically decreasing function of the shear Reynolds 



number, which is analogous to the reduction of rric by 
shear that we observed. Transverse rolls also travel un- 
der unidirectional plane Couette shear, again analogous 
to our travelling vortex state. 

Whereas unbounded RBC with plane Couette shear 
cannot be realized experimentally, RBC has been stud- 
ied experimentally and theoretically in narrow slots with 
open through- flows. [^ ^4|j45|] Quasi-2D transverse rolls 
can be stabilized by wall effects in slots. The through- 
flow consists of a weak Poiseuille flow with a very small 
Reynolds number. Its effects on RBC are well under- 
stood. In brief, the onset of convection is again sup- 
pressed, but the first instability is convective (i.e. it 
grows only downstream of a localized perturbation), 
rather than absolute. The resulting convection pattern 
drifts in the direction of the through flow. It is inter- 
esting that the clear distinction between convective and 
absolute instability is blurred in annular electroconvec- 
tion with shear, in which the 'through' flow loops back 
on itself. The annular geometry is naturally closed. The 
linear stability analysis of Ref. only treated absolute 
instablity. In principle, the system might still be consid- 
ered to be convectively unstable if localized perturbations 
grew only as they travelled azimuthally. 

Taylor Vortex flow (TVF) [Q is an extensively studied 
pattern-forming instability with geometric similarities to 
annular electroconvection. However, the instability lead- 
ing to TVF depends crucially on axial disturbances to 3D 
Couette flow. Purely 2D circular Couette flow is, in fact, 
linearly stable. ^6 47 1 TVF is the result of an instability 
of a 3D shear flow, while what we have studied here is the 
effect of a shear flow on the electroconvective instability. 

Agrait and Castellanos theoretically studied the effect 
of a 3D Couette shear on an electrohydrodynamic insta- 
bility in TVF geometry. Their electrohydrodynamic 
system consisted of a nearly insulating fluid confined be- 
tween metallic electrodes. Charge injection, a process by 
which charge carriers are created at the electrodes, occurs 
when strong electric fields are applied. The interaction of 
this volume charge density with the applied electric field 
leads to electroconvection instabilities. |lj,|l^,|l^ Agrait 
and Castellanos considered electroconvection due to a ra- 
dial field with charge injection on either cylinder. Both 
cylinders were permitted to rotate to produce a general 
Couette shear. They found that shearing enhanced the 
instability, leading to a 3D flow that resembles TVF. This 
is in direct contrast to the stabilizing effects we observed 
in 2D. 

Although shear and rigid rotation are conceptually 
quite different, it is interesting to compare their effects on 
instabilities. Again, we flnd crucial distinctions between 
2D and 3D systems. 

The added stability in sheared annular electroconvec- 
tion is a consequence of the shear and not of rotation. 
Under rigid rotation, where the inner and outer elec- 
trodes are co-rotating, one can transform to rotating co- 
ordinates in which the electrodes are stationary. This 
transformation introduces a Coriolis term —2Uz x v = 
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— 2r2V'0 in the Navier-Stokes equation which may be ab- 
sorbed into the pressure gradient term VP and eUmi- 
nated. Thus, in a purely 2D system, rigid rotation 

and the non-rotating, unsheared case have identical sta- 
bility. It also follows, since the transformation is gen- 
eral and the unsheared bifurcation is stationary, that the 
resulting nonlinear vortex pattern above onset must be 
stationary in the co-rotating frame. 

This lack of dependence on rigid rotation may be con- 
trasted with a large class of 3D and quasi-2D rotating 
Rayleigh-Benard systems P,p8|-^ , where rotation pro- 
duces added stability but the absence of strictly 2D flow 
results in a time-dependent (precessing) convection pat- 
tern in the co-rotating frame. The precession seen in 
these systems is analogous to the travelling of the sheared 
patterns that we observed only to the extent that rota- 
tion and shear break the same symmetry, which allows 
travelling patterns. Symmetry aside, the physical ori- 
gins of added stability and precession are fundamentally 
different between the 2D sheared and 3D rotating cases. 

In general, 3D rotating convecting systems may also 
support strictly 2D 'Taylor column' solutions which 
do not precess in the rotating frame and whose onset 
occurs at the same critical Rayleigh number as in the 
absence of rotation. [^||30|] The system most analogous 
to ours is the interesting but experimentally unrealiz- 
able situation of 2D RBC in a rotating annular geometry 
with purely radial gravity and heating |^o| . Theoretical 
studies found columnar solutions which are very closely 
analous to our vortices. One might hope to approach this 
limit in experiments using radially temperature gradients 
imposed between rapidly rotating concentric cylinders. 
Purely columnar solutions have not observed in any ro- 
tating RBC experiment because the boundary conditions 
at the top and bottom of the cylinder must be stress free 
iQ,^, a requirement that cannot be attained in terres- 
trial RBC experiments. In contrast, two-dimensionality, 
stress free surface conditions and radial driving forces 
all arise naturally in the electroconvection of an annular 
suspended film. Thus, the vortices that occur in annular 
electroconvection without shear are accurate 2D analogs 
of 'Taylor columns'. 

The similarity between our system and many other 
better-studied systems and models suggests that many 
linear and nonlinear techniques developed for other prob- 
lems can be fruitfully brought to bear on annular electro- 
convection. In addition, the system may allow the study 
of bifurcation scenarios that are not experimentally real- 
izable in other similar systems. 



VI. CONCLUSION 

In this paper, we have reported a wide ranging exper- 
imental study of the primary bifurcation to electrocon- 
vection in sheared two-dimensional annular films. Our 



principal experimental probe was the excess current car- 
ried by the film due to convection, which can be directly 
related to the amplitude of the convective flow. In all, we 
examined annuli with six different radius ratios a, with 
0.33 < a < 0.80. For these, the Reynolds number of 
the applied Couette shear varied between < TZe < 3. 
The explored range for the dimensionless parameter V, 
which varies with film thickness and conductivity, was 
1 < P < 150. 

We compared this data, in the first instance, with the 
predictions of linear theory |3j. The data for the criti- 
cal voltage for the onset of convection without shear 
could be combined for all six radius ratios and was shown 
to obey the scaling predicted by linear theory. This data 
could then be used in a slightly model-dependent way to 
fix the only remaining unknown material parameter, the 
fluid viscosity. We could then compare the experimen- 
tally measured suppression of the onset by shear to linear 
theory in an essentially pararmeter-free way. The agree- 
ment was satisfactory for a wide range of a, V, and TZe. 
Using a simple visualization scheme, we also confirmed 
that the azimuthal mode number m near onset was close 
to the most unstable critical value rUc predicted by linear 
theory. 

Nonlinear fits to current-voltage data above the on- 
set of convection were then used to infer the real part 
of the amplitude of convection. Under shear, the ampli- 
tude can be shown to obey a complex Ginzburg-Landau 
equation; its real part was fit to the time-independent 
steady state version of this equation in order to extract 
the coefficient g of the cubic nonlinearity. When shear 
is absent, the primary bifurcation is a pitchfork while 
with shear, the bifurcation is a pitchfork Hopf. We de- 
termined the functional dependence of g as a, TZe and T' 
were varied. When shear is absent, it was found that for 
a = 0.47, 0.56, 0.60, 0.64, 0.80 and P ^ 13 that g was in- 
dependent of "P. In addition, in this regime g > 0, hence 
the bifurcation was supercritical. For a — 0.33, g was 
found to be an increasing function of "P for 2 < P < 8. 
More importantly, it was found that the bifurcation was 
subcritical((? < 0) for P ~ 5 and supercritical(g > 0) for 

P ^ b. In overall trend, we found that g is an increasing 
function of a for zero shear. The largest a — 0.80 was 
sufficiently close to the limiting case of a ^ 1 that we 
could extrapolate the measured values of g for compar- 
ison to the value of g calculated from weakly nonlinear 
theory for the corresponding 'plate' electrode geometry. 
This quantitative comparison gave reasonable agreement. 

When shear was applied, it had a strong effect on 
the subcriticality of the primary bifurcation. Measure- 
ments of 5 as a function of the shear Reynolds number 
TZe revealed that for a = 0.47, 0.56, 0.60, 0.64, 0.80 there 
existed a tricritical Reynolds number TZer below which 
g > and the onset is a supercritical Hopf bifurcation. 
For TZe > TZer, 5 < and the onset became a subcrit- 
ical Hopf bifurcation with the degree of subcriticality a 
nontrivial function of TZe. Hence, the Reynolds number 
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forms an interesting second control parameter in this sys- 
tem which can be used to vary the nature of the primary 
bifurcation. 

The data presented in this paper represents only 
a sparse sampling of the full 3D parameter space of 
g{a, V,TZe). Within this space, there presumably exist 
continuous 2D loci on which g ~ and the primary bi- 
furcation is tricritical. On other loci, the hysteresis Se is 
a local maximum. As a function of a, linear theory |^] 
shows that there exist special values of a at which two 
aziniuthal mode numbers m and m+l are simultaneously 
linearly unstable. All of the g results we have described 
pertain to the primary bifurcation; above onset we also 
observed numerous secondary bifurcations which take the 
form of hysteretic jumps in m. The location of these sec- 
ondary bifurcations is a function of TZe. These consider- 
ations suggest that more-or-less discontinuous jumps in 
the behavior of g may occur at parameter values across 
which the azimuthal mode number m of the nonlinear 
pattern changes discontinuously. The rich phenomenol- 
ogy of even the primary bifurcation presents a significant 
challenge to the weakly nonlinear theory of this instabil- 
ity. 
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TABLE I. Experimental measurements of the marginally 
stable mode number for zero shear, mc-Theoretical values are 
from Ref. [3]. 
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TABLE II. Experimental measurements of the coefficient 
of the cubic nonlinearity, g without shear. 
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TABLE III. Experimental measurements of the tricritical 
Reynolds number TZbt at which g = 0. 



radius ratio 
a 


Minimum 
9 


Reynolds number 


r 


Maximum 
Tie 


0.47 


-3.68 ± 0.19 


0.83 ± 0.18 


15.3 


2.59 


0.56 


-5.15 ± 1.04 


0.11 ± 0.05 


63.3 


0.22 


0.60 


-1.74 ± 0.04 


0.05 ± 0.02 


32.1 


0.13 


0.64 


-4.34 ± 0.79 


0.23 ± 0.02 


53.4 


0.25 


0.80 


-9.17 ± 0.56 


0.04 ± 0.01 


12.0 


0.10 



TABLE IV. Experimental measurements of the minimum 
value of g and the corresponding TZcmin and P values. 
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